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Abstract 

In a recent paper [Nature 412, 712 (2001)], Zurek has argued that 
(1) time evolution typically causes chaotic quantum systems to generate 
structure that varies on the scale of phase-space volume elements of size 
(h 2 /A) d , where A is a classical action characteristic of the state and d is 
the number of degrees of freedom, and that (2) this structure implies that 
a small change in a phase-space coordinate X by an amount SX ~ HX/A 
generically results in an orthogonal state. While we agree with (1), we 
argue that (2) is not correct if the number of degrees of freedom is small. 
Our arguments are based on the Berry-Voros ansatz for the structure 
of energy eigenstates in chaotic systems. We find, however, that (2) 
becomes vàlid if the number of degrees of freedom is large. This implies 
that many-body environments may be crucial for the phenomenon of 
quantum decoherence. 
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Zurek |J has recently pointed out that, contrary to prior expectations @, a quantum 
system with d degrees of freedom will generically develop structure on the scale of phase- 
space volume elements of size (h 2 /A) d , where A is a typical classical action associated with 
the state. This is diagnosed by considering the Wigner function for the quantum state 

/d d s 
j—^ exp(ip-s/h)íp*(x + js)ífj(x - ±s), (1) 

where ip(x) = (x\ip) is the position-space wave function. (We have assumed a pure state; the 
gener alization to mixed states is given in [1].) By considering some simple superpositions 
of gaussian wave packets (coherent states), Zurek shows that there can be fluctuations of 
W(x,p) in the position x on scales A Ximin ~ h/P (where P is a typical classical value of the 
momentum p), and in p on scales X p , 

min ~ h/L (where L is a typical classical value of x^j. 
The corresponding phase-space volume is (A^minÀpmin)^ ~ (h 2 /A) d , where A ~ PL. 

Zurek goes on to consider the possibility that this sub-Planck structure may have im- 
portant physical significance if the system in question serves as an environment for a qubit 
(or other two- or few-state system). In particular, he argues that if the two possible states 
of the qubit cause slightly different displacements in the phase space of the environment, 
then these two possible environmental states will be nearly orthogonal. This is then tanta- 
mount to decoherence of the quantum state of the qubit by the environment. To be specific, 
consider the action of a unitary displacement operator 

D(5p,5x) = exp[i(5p-x + Sx-p)/h\, (2) 

where x and p are the position and momentum operators, on a state \ip) of the environment. 
The claim of [1] is that 

(ip\D(5p,8x)\ip) ~0 for 6x>h/P or 5p>h/L. (3) 

We will argue, however, that this is not correct in general for systems with few degrees of 
freedom. On the other hand, we find that eq. (||]) is generically correct if the system (which 
is playing the role of the environment) has many degrees of freedom. 
Much of the analysis of [1] is based on the formula 

\(ilj\D(óp,5x)\tlj}\ 2 = (2irh) d J d d xd d pW(x,p)W(x+5x jP +5p). (4) 

However, a somewhat simpler and equally vàlid formula is 

(i/)\D(6p,6x)\i/)) = J d d xd d pé {5p - x+5x - p)/h W{x,p). (5) 

Here we have made use of one of the bàsic properties of the Wigner function: expectation 
vàlues of symmetrically ordered operators can be written as simple integrals. Note that all 
x's and p's in eq. (0) are now classical integration variables. Also, in the special case 5x = 0, 
eq. (HD reduces to eq. (16) of [1]. 

From eq. (|5p, we see immediately that we do not need to have any information about 
the small-scale structure in W(x,p) in order to draw conclusions about the behavior of 
(D(5p,5x)) for small 5p and Sx. This follows simply from the mathematics of the Fourier 
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transform, which relates small-scale structure in the original variables (x and p) to large- 
scale structure in the transformed variables (óp and óx). In fact, according to eq. ([|), the 
implication of the small-scale structure in W(x,p) is that (D(5x,5p)) will be generically 
nearly zero only for óx > ^/A Pimin ~ L or óp > h/\ x ,mm ~ P. That is, to get a generically 
orthogonal state, the shift in coordinates or momenta must be comparable to the classical 
scales that characterize the original state. 

On the other hand, what happens to (D(óp,óx)) at small óp or óx is governed in turn 
by the behavior of W(x,p) for large x or p. If we suppose that the corresponding state has 
spread out over the available phase space (which is characteristic of chaotic systems, and is 
explicitly assumed in [1]), then we can treat the large-scale structure of the Wigner function 
via the Berry-Voros ansatz f||-|5| 

W(x,p) ~ ó(H(x,p) — E), (6) 

where ó(E) is the Dirac delta function, H(x,p) is the hamiltonian, and E is the energy 
(assumed to be classically sharp). This ansatz simply corresponds to the ergodic distribution 
that covers the constant-energy surface uniformly, according to the Liouville probability 
measure. It is devoid of the small-scale structure elucidated in [1], but correctly describes 
the large-scale structure. When eq. (|5|) is used in eq. (Q), the restriction of the range of 
integration to a thin shell of phase space will produce ringing in (D(óp,óx)). The first 
zeroes will occur at óx ~ h/P and óp ~ h/L; these are the scales where, according to [1], 
(D(óp,óx)) should become small. This is correct, but then the ringing effect implies that 
there are recurrent màxima and minima in (D(óp,óx)), with only a power-law fall-off of 
the peak heights. (The exponents in these power laws depend on the details of the system; 
more on this below.) This is contrary to the supposition of [1] that any recurrences will be 
strongly suppressed in general (even though they occur in the example of four superposed 
gaussians that is worked out explicitly). Finally, for óx ~ L and óp ~ P, the small-scale 
structure of W(x,p) becomes important, and (D(óp, óx)) goes rapidly to zero. 

To illustrate this, we consider a paradigmàtic chaotic system, a two-dimensional billiard 
such as the stadium ||. In this case H(x,p) = ^p 2 when x is inside, and H(x,p) = oo 
when x is outside. If (for pedagogical simplicity) we take the billiard to be nearly circular 
with radius L (e.g, a stadium with a short straight segment 0, or a circle with rough walls 
||), and the bouncing ball to have energy E = ^P 2 , then from eqs. (|5|) and (^) we get 

(D(óp, óx)) = J (P\óx\/h) Jl ^j£\ ( 7 ) 

where J n (£) is an ordinary Bessel function. Here the oscillations and slow decay of the Bessel 
functions correspond to the ringing phenomenon that we described above. Also, because 
the Berry-Voros ansatz does not have the correct small-scale structure, we do not see the 
cutoffs at óp ~ P and óx ~ L that would otherwise be there. 

We can also go a little further than this for a chaotic system. Some time ago, Berry 
H conjectured that, in such a system, the energy eigenfunctions (in position space) have 
the statistical properties of a gaussian-random superposition of plane waves, all with the 
wavenumber k given locally by classical energy conservation and the fundamental quantum 
relation k = P/h, but with wave- vectors pointing in random directions. Since then, this 
conjecture has been supported by a number of numerical experiments [pHT^I, and it underlies 
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the successful characterization of the distribution of conductance peak heights in quantum 
dots in the Coulomb-blockade regime [13-fLql. If we use an overbar to denote the statistical 



random-wave average, then Berry's conjecture can be neatly summarized via 



W(x,p) ~6{H{x,p)-E), (8) 

plus the statistical statement that the wave-function ip(x) is a gaussian random variable. 
This in turn can be shown to imply [fT^| 



W(x,p) 2 = W(x, P y + 0{h d ~ l ) (9) 



for d > 2. In this formula, W(x,p) itself is to be regarded as 0(h°). Eq. @ means that 
each particular energy eigenfunction will have only small deviations from eq. (§) on large 
scales. Furthermore, eqs. and (|9|) also apply for genèric time-dependent superpositions 
of energy eigenstates, 

m)) = Y,c«e- lEat/h \a), (10) 



and for genèric mixed states, 



p(t) = E^^ { " " W >)(Pl (11) 

a/3 

provided we assume that different eigenfunctions are statistically independent superpositions 
of random waves. (Also, we assume that the quantum energy uncertainty is classically small.) 
Thus, in order to construct a state that is, for example, localized in a particular small region 
of phase space, we must carefully choose the magnitudes and phases of the c a 's in eq. ( |Ï0f ) 
or the c aj g's in eq. (fïl"|). Once we have done this, we expect that subsequent time evolution 
will eventually spread the packet over the corresponding energy surface, its distribution on 
large scales governed by the classical ergodic distribution. In quantum mechanics, this time 
evolution is given entirely by the phase changes in eq. ([Ï0|), and so these must be capable of 
producing this spreading. This is guaranteed by eqs. (^) and (Q). In a many-body system, 
this same effect corresponds to the approach to thermal equilibrium fj~6|--|T8|. 



We have shown that the small-scale structure in phase space discovered in [1] does not 
generically result in small vàlues of the overlap between a displaced state and the original 
unless the displacements are classically large. However, this conclusion does not take into 
account new effects that can arise when d, the number of degrees of freedom, is large. 
Essentially, what happens for d 3> 1 is that the power-law exponent for the fall-off of the 
recurrent peaks in (D(5p,5x)) becomes large with d, suppressing all the recurrences, and 
restoring the conclusions of [1]. We now turn to our exposition of this phenomenon. 

For defmiteness, we will consider a dilute gas of N hard spheres in a three-dimensional 
box of volume L 3 . This system can also be viewed as a single bouncing ball in a complicated 
3iV-dimensional chaotic billiard, and so we will assume that Berry's random-wave conjecture 



holds for it at sufficiently high total energy E |Ï6|. We take E = ^P 2 , so that P is the 



typical momentum of a single gas particle. The appropriate generalization of eq. (^) is then 
(D(6p Sx)) - UNV*P\óx\/h) fí sin(L5 Pl /2h) 

{U[dp ' dX)) (N^P\Sx\/2hyv\ l\ LS Pi /2h ' [ ' 
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where v = (3N — 2)/2. To perforin the position-space integral that results in the múlti- 
ple product, we ignored the spatial regions that are excluded by the requirement that the 
hard-sphere centers not come closer than the hard-sphere diameter; this should not be an 
important effect for a dilute gas. In the large-iV limit, the Sx dependent factor in eq. (|Ï2|) 
becomes exp(— P 2 \Sx\ 2 /6h ). If we assume that the Sp^s are all of comparable magnitude, 
and that each is less than h/L, then the Sp dependent factor becomes exp(— L 2 \5p\ 2 /24h 2 ). 
Thus we have 

(D(5p, 5x)) = exp{-P 2 \5x\ 2 /Qh 2 ) exp{-L 2 \5p\ 2 /24h 2 ). (13) 

From this formula we see that (D(5p,5x)) goes rapidly to zero for either \5x\ > h/P or 
I^pI ftfi'/L, which is now in accord with [1]. 

Before concluding, we briefly return to systems with few degrees of freedom, and consider 
the reconciliation of eqs. (^) and (Bl). Assuming that the state spreads out over the available 
phase space, we have seen that eq. © implies that (D(5p, Sx)) rings, with recurrent zeroes at 
múltiples of Sx ~ h/P and óp ~ h/L. How can this behavior be reproduced by eq. (^)? It is 
clear that this is possible if and only if W(x, p) in fact has the small-scale structure identified 
in [1]. Thus there is an intimate relation between the large- and small-scale structure of the 
Wigner function. 

To summarize, we find that small displacements 5x > h/P in position and/or 5p > h/ L of 
the quantum state of a many-body chaotic system generically produce a nearly orthogonal 
state, provided the system has many degrees of freedom. For a chaotic system with few 
degrees of freedom, we find that the overlap between the displaced and the original state 
has a series of zeroes spaced by 5x ~ h/P and 5P ~ h/L, separated by broad peaks whose 
heights fall off as a power of 5x or 5p, with a final exponential cutoff at 5x ~ L and óp ~ P. 
These results imply that a many-body system, thought of as an environment, will (not 
unexpectedly) be more effective in causing a simple system (such as a qubit) to undergo 
decoherence. 
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